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A  simulation  of  alternating-current  thin-film  electroluminescent  device  operation  with  positive  space 
charge  present  in  the  phosphor  layer  of  the  device  is  presented.  The  simulation  is  based  on  modeling 
the  space-charge  distribution  using  a  single-sheet  charge  model.  The  simulation  is  performed  for 
two  cases  of  space-charge  creation:  by  impact  ionization  of  deep  levels  in  the  phosphor  or  by  field 
emission  from  traps  in  the  phosphor.  Results  of  the  simulation  show  that  space-charge  creation  by 
either  mechanism  is  capable  of  causing  overshoot  in  both  capacitance-voltage  and  internal 
charge-phosphor  field  ( Q-Fp )  plots.  ©  1995  American  Institute  of  Physics. 


I.  INTRODUCTION 

A  common  simplifying  assumption  in  many  previously 
developed  alternating-current  thin-film  electroluminescent 
(ACTFEL)  device  models  is  that  the  entire  phosphor  layer  of 
the  device  is  a  constant  field  region.1-5  The  assumption  of 
constant  phosphor  field  implies  that  bulk  space  charge  does 
not  exist  in  the  phosphor  layer.  In  modeling  ZnS:Mn  devices 
grown  by  evaporation,  the  results  obtained  under  the  as¬ 
sumption  of  a  constant  phosphor  field  match  the  measured 
results  very  well.6,7  However,  with  the  advent  of  alternative 
phosphor  materials  and  methods  of  fabrication,  certain  ef¬ 
fects  have  been  experimentally  observed  that  cannot  be  de¬ 
scribed  using  a  constant  phosphor  field  model. 

The  effects  commonly  observed  in  experimental  mea¬ 
surements  of  ACTFEL  devices  that  are  not  consistent  with 
the  constant  phosphor  field  model  are  brightness-voltage 
( B-V )  hysteresis,2,8"10  negative  differential  resistance,11,12 
capacitance-voltage  (C-V)  overshoot,13,14  and  internal 
charge-phosphor  field  ( Q-Fp )  overshoot.13-15  Other  re¬ 
searchers  have  developed  numerical  models  to  account  for 
the  existence  of  space  charge  in  the  phosphor  layer  of  an 
ACTFEL  device,  primarily  focusing  on  the  issues  of  B-V 
hysteresis9,16-18  and  negative  differential  resistance.11,12  In 
contrast,  the  work  presented  herein  is  concerned  with  mod¬ 
eling  both  C-V  and  Q-Fp  overshoot  using  a  single-sheet 
charge  model7  to  account  for  space  charge  in  the  phosphor. 
In  particular,  effects  seen  in  ZnS:Mn  devices  grown  by 
atomic  layer  epitaxy  (ALE)  are  considered.  The  information 
gained  from  simulation  of  the  single-sheet  charge  model  is 
then  used  to  aid  in  the  interpretation  of  laboratory  measure¬ 
ments  of  ACTFEL  devices  exhibiting  space-charge  related 
effects. 

The  two  distinguishing  characteristics  of  the  work  pre¬ 
sented  herein  is  the  utilization  of  a  single-sheet  charge  model 
and  the  application  of  this  modeling  towards  obtaining  an 
understanding  of  C-V  and  Q-Fp  overshoot.  The  single-sheet 
charge  ACTFEL  device  model  is  chosen  over  distributed 
models11,12,19  for  two  primary  reasons.  First,  the  lumping  of 
the  space-charge  distribution  into  a  single  sheet  greatly  sim¬ 
plifies  the  problems  associated  with  the  presence  of  space 
charge  in  the  phosphor  layer.  Second,  the  highly  polarity- 


dependent  electrical  characteristics  of  ALE  ZnS:Mn  devices 
being  modeled  suggest  that  the  space-charge  distribution  is 
concentrated  near  the  edges  of  the  phosphor  layer,  rather 
than  evenly  distributed  across  the  width  of  the  phosphor 
layer. 

It  should  also  be  noted  that  most  of  the  previous  research 
performed  on  modeling  space-charge-related  effects  of  ACT¬ 
FEL  devices  has  assumed  the  method  of  space-charge  cre¬ 
ation  is  the  trapping  of  holes  generated  by  band-to-band  im¬ 
pact  ionization  processes.9,12,17,19  However,  Yang  et  al 
provide  experimental  evidence  that  hot-electron  impact  ion¬ 
ization  of  deep  levels  is  more  likely  to  cause  hysteresis  than 
hole  trapping.10  Therefore,  the  focus  of  the  modeling  in  this 
paper  is  towards  modeling  space-charge  creation  as  due  to 
either  impact  ionization  of  deep  levels  or  field  emission  of 
electrons  from  bulk  traps. 

II.  SINGLE-SHEET  CHARGE  MODEL 
A.  Single-sheet  charge  model  quasistatics 

The  single-sheet  charge  model  for  an  ACTFEL  device 
models  a  space-charge  distribution  in  the  phosphor  layer  as  a 
sheet  of  charge  at  a  single,  arbitrary  location  within  the  phos¬ 
phor  layer.7  The  general  form  of  the  single-sheet  charge 
model  is  depicted  in  Fig.  1.  In  Fig.  1,  ix  and  i2  refer  to  the 
ACTFEL  insulator  layers,  the  subscript  “sc”  refers  to  the 
single-sheet  charge  layer,  and  px  and  p2  refer  to  the  two 
different  regions  of  the  phosphor  within  which  the  electric 
field  may  be  different  because  of  the  existence  of  the  space- 
charge  layer.  The  labels  /,  d ,  and  q  refer,  respectively,  to  the 
electric  field,  charge  density,  and  thickness  of  the  subscripted 
layers. 

To  develop  an  analytical  description  for  the  single-sheet 
charge  model  to  be  simulated,  various  basic  principles  are 
applied  to  derive  a  set  of  equations  that  describe  the  quasis¬ 
tatics  of  the  model.  This  process  begins  by  application  of 
Kirchoff’s  voltage  law  to  the  structure  of  Fig.  1  to  yield 

dnfn(t)  +  dsfpi(t)  +  {dp-ds)fp2(t)  +  di2fi2(t)=-vg(t). 

(1) 

Two  things  should  be  noted  from  the  above  equation  and 
from  Fig.  1.  First,  the  sign  convention  of  the  electric  field 
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FIG.  1.  Single- sheet  charge  model  of  an  ACTFEL  device. 


terms  in  Eq.  (1)  is  opposite  the  usual  textbook  convention  for 
electric  fields.  This  convention  is  adopted  to  stay  consistent 
with  the  formulation  presented  by  Bringuier3  and  other  work 
based  on  the  use  of  this  convention.13-15  Second,  the  external 
series  capacitor  Cs  is  neglected  in  formulating  Eq.  (1).  This 
is  a  valid  assumption  because  in  laboratory  measurements  Cs 
is  chosen  to  be  large  compared  to  the  capacitance  presented 
by  the  ACTFEL  device.  Next,  Gauss’  Law  is  applied  to  each 
of  the  charge  interfaces  to  yield  the  following  three  equa¬ 
tions: 

cndnfn(t)-cpdpfpl(t)=-q1(t),  (2) 

cpdp[fpi(t)-fp2(t)]=-qsc(t),  (3) 

cpdpfpl{t)-ci2di7fn{t)=-q2(t),  (4) 

where  the  c  terms  represent  the  capacitance  of  the  sub¬ 
scripted  layer.  Then,  it  is  necessary  to  add  an  external  equa¬ 

tion  that  relates  the  charge  on  the  initially  cathodic  electrode 
to  the  charge  measured  on  a  capacitor  in  series  with  the 
ACTFEL  device.  This  equation  represents  the  bridge  be¬ 
tween  experimentally  measurable  information  and  the  inter¬ 
nal  device  physics.  The  external  equation, 

ciidnfn(t)=  (5) 

remains  unchanged  from  the  model  presented  by  Bringuier 
that  assumes  no  space-charge  creation  in  the  device.3  Finally, 
invoking  the  necessity  for  charge  balance  inside  the  device,  it 
is  known  that 

q\(t)  +  q2{t)  +  q!iC(t)  =  0.  (6) 

The  above  six  equations  serve  to  completely  describe  all 
charge  densities  and  fields  in  the  single-sheet  charge  model. 
The  scope  of  their  utility  will  become  evident  throughout  the 
development  of  the  dynamic  single-sheet  charge  model. 

The  first  steps  in  the  development  of  the  dynamic  single¬ 
sheet  charge  model  are  to  perform  several  manipulations  of 
Eqs.  (1)  through  (6).  First,  adding  Eqs.  (2)  through  (4)  and 
substituting  Eq.  (6)  into  the  sum  yields  the  following  result: 

ci\di\fn(t)  =  ci2di2fi *2(0*  0) 

Since  cndn  and  ci2di2  represent  the  dielectric  permittivities 
of  each  respective  insulator  layer,  it  follows  that  the  insulator 
fields  divide  according  to  the  ratio  of  the  dielectric  constants 
of  the  respective  insulating  materials.  Therefore,  for  ACT¬ 
FEL  devices  with  insulating  layers  of  the  same  material  and 


thickness,  fn{t)—fi2{t).  Throughout  the  remainder  of  this 
article  it  is  assumed  that  the  device  being  modeled  is  fabri¬ 
cated  such  that  both  insulator  layers  are  identical. 

The  next  step  towards  obtaining  the  dynamic  model  for 
the  evolution  of  the  phosphor  fields  is  to  find  expressions  for 
both  fpl  and  fp2.  When  Eqs.  (2),  (4),  (5),  and  (7)  are  solved 
for  fpl  and  fp2,  respectively,  the  expressions 

[<?!«-<?*(*)],  (8) 

cpap 

fP2(t)=--—rlq2(t)  +  <le(t)]  (9) 

Lpap 

are  obtained.  The  significance  of  Eqs.  (8)  and  (9)  lies  in  the 
fact  that  expressions  are  obtained  for  fp\(t)  and  fp2(t)  that 
depend  exclusively  on  the  values  of  the  three  charge  vari¬ 
ables,  qx(t)9  q2(t),  and  qe(t). 


B.  The  dynamics  of  the  single-sheet  charge  model 


To  obtain  the  basis  of  the  dynamic  model,  the  expres¬ 
sions  for  fpi  and  fp2  in  Eqs.  (8)  and  (9)  are  differentiated 
with  respect  to  time.  The  resultant  equations  are  the  follow¬ 
ing  pair  of  coupled,  first-order  differential  equations: 


<?/,!„  1  /<?<?! 

dt  Cpdp  \  dt  dt  /’ 


(10) 


df  p2  _  __  1  /  ^2  +  dcte\ 

dt  Cpdp  \  dt  dt  I 


(ID 


Next,  when  Eqs.  (1)  through  (6)  are  self-consistently  solved 
for  qe  in  terms  of  qx ,  qsc,  and  vg  ,  and  the  result  is  differen¬ 
tiated  with  respect  to  time, 


d<le 

dt 


Ci 

'*>+( 

■4) 

dpl 

dqx 

dvs 

^  _ £L 

Ct+Cp 

dt  \ 

dt 

i  C  f 

t  dt 

(12) 


is  obtained.  The  expression  for  dqjdt  given  by  Eq.  (12)  and 
the  charge  neutrality  condition  given  by  Eq.  (6)  are  then 
substituted  into  Eqs.  (10)  and  (11)  to  give  expressions  for 
dfp\f  dt  and  dfp2Idt  in  terms  of  the  time  derivatives  of  the 
charge  on  each  interface.  The  equations 


dfpi .  1 

dqi  1  ct  \ 

fx_ 

*l) 

dq*c 

dt  cpdp 

{  dt  \c<  +  cp/ 

H 

l1 

dp) 

dt 

+  ct 


dVg 

dt 


(13) 


dfp2  _  1 

'  dq2  (  C,-  \ 

dq2 

lds\ 

<^7  SC 

dt  Cpdp 

[  dt  \  Cj~h  CpJ 

17  +  l 

[ dpi 

dt 

+  cf 


to*] 

dt  J 


(14) 


result.  The  above  pair  of  differential  equations  is  complete  in 
that  it  describes  the  dynamics  of  the  phosphor  fields  in  terms 
of  the  three  internal  current  density  terms,  dqfdt ,  dq2ldt, 
and  dqjdt  along  with  the  external  voltage  slew  rate  term, 
dvg!dt .  To  complete  the  dynamic  model,  it  is  necessary  to 
develop  expressions  for  dqfdt ,  dq2ldt ,  and  dqjdt .  To  de- 
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FIG.  2.  Energy  band  diagram  depicting  the  four  space-charge  current 
components. 


velop  expressions  for  these  three  terms,  it  is  first  necessary  to 
examine  the  situation  at  each  charge  interface  more  closely. 

The  first  interface  charge  to  consider  is  the  sheet  of 
space  charge  located  at  some  arbitrary  location  in  the  phos¬ 
phor  bulk.  The  space-charge  layer  is  depicted  in  Fig.  2  in  an 
energy  band  diagram  of  the  phosphor  layer  of  the  ACTFEL 
device.  As  seen  in  Fig.  2,  the  space-charge  layer  can  emit 
electrons  to  either  interface  1  or  interface  2,  depending  on 
the  field  polarity  in  the  phosphor  layer.  Also,  the  space- 
charge  layer  can  trap  electrons  emitted  from  either  interface 
1  or  interface  2.  Each  of  the  electron  emission  or  reception 
events  described  above  now  must  be  converted  to  equation 
form  to  describe  the  time  derivative  of  qsc.  To  accomplish 
this,  usual  current  conventions  are  employed;  i.e.,  an  electron 
emitting  interface  corresponds  to  a  positive  current  into  the 
interface  in  question,  whereas  a  receiving  interface  corre¬ 
sponds  to  a  negative  current  into  the  interface  in  question. 
Next,  the  field  dependencies  for  electron  emission  need  to  be 
deduced.  It  is  seen  in  Fig.  2  that  both  the  emission  of  elec¬ 
trons  from  interface  1  to  the  space-charge  layer  and  the  emis¬ 
sion  of  electrons  from  the  space-charge  layer  to  interface  1 
are  controlled  by  the  field  magnitude  and  polarity  in  phos¬ 
phor  region  1  (pj).  Similarly,  electron  emission  to  and  from 
interface  2  is  controlled  by  the  field  magnitude  and  polarity 
in  phosphor  region  2  (p2).  Finally,  from  the  above  consider¬ 
ations,  the  equation  that  is  sought  emerges, 

= JJJp Aeak(fp i )  - ,  (is) 

where  J  terms  represent  current  densities  from  the  sub¬ 


scripted  emitting  interface.  It  should  be  noted  that  the  super¬ 
scripts  on  the  latter  two  terms  in  Eq.  (15)  represent  the  cur¬ 
rent  emitted  by  the  subscripted  semiconductor-insulator 
interface  that  annihilates  positive  space  charge.  This  current 
is  denoted  /leak  because  it  is  standard  to  call  charge  emitted 
during  the  interpulse  interval  leakage  charge,13"15  and  space 
charge  is  mostly  annihilated  during  the  low-field  portion  of 
the  interpulse  interval. 

Continuing  the  same  approach  outlined  above,  Fig.  3  can 
be  employed  to  produce  a  similar  expression  for  dqxldt. 


FIG.  3.  Energy  band  diagram  depicting  the  current  components  of  inter¬ 
face  1. 

First,  with  the  correct  field  polarity  in  the  phosphor,  electrons 
are  emitted  from  interface  1.  Then,  with  the  opposite  field 
polarity,  electrons  are  emitted  by  the  space-charge  layer  and 
received  by  interface  1.  Also,  concomitant  with  the  space- 
charge  emission,  electrons  are  emitted  from  interface  2. 
However,  some  of  these  electrons  are  captured  by  the  space- 
charge  layer  and  never  reach  interface  1.  Therefore,  the  su¬ 
perscript  term  from  the  development  of  the  expression  for 
dqjdt  is  the  term  that  accounts  for  this.  The  total  current 
reaching  interface  1  is  then  the  difference  between  J2  and 
yi,eak.  This  is  shown  pictorially  in  Fig.  3.  Finally,  the  expres¬ 
sion  for  dq xIdt  is 

^ =Jl(fpl) -Ufpl) -Jscifpl)  +  -/2ak(/p2).  (16) 

where  the  J  terms  again  represent  emission  from  the  sub¬ 
scripted  interface.  The  Jsc  and  are  the  same  terms  used 
in  the  development  of  the  expression  for  dqjdt. 

To  complete  the  set  of  expressions  for  substitution  into 
Eqs,  (13)  and  (14),  a  similar  expression  for  dq2ldt  must  be 
developed  next.  The  figure  of  interest  for  electron  emission 
and  reception  by  interface  2  is  Fig.  4.  The  development  of 
this  expression  exactly  parallels  the  development  of  the  ex¬ 
pression  for  dqxldt ,  so  the  details  will  be  omitted.  The  de¬ 
sired  expression  is 

~  =  hifpi)  -  Jx  (fP  i )  -  JJfpi)  +  •  (17) 

Equations  (15)  through  (17)  generated  above  comprise  the 
complete  set  of  equations  necessary  for  substitution  into  Eqs. 
(13)  and  (14). 

III.  ELECTRON  EMISSION 

To  complete  the  simulation  model,  expressions  must  be 
developed  that  describe  each  of  the  current  density  terms  in 
Eqs.  (15)  through  (17).  The  approach  to  accomplishment  of 
this  goal  is  to  first  examine  electron  emission  from  the 
insulator-phosphor  interfaces.  Then,  once  this  is  established, 
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FIG.  4.  Energy  band  diagram  depicting  the  current  components  of  inter¬ 
face  2. 


space-charge  creation  is  examined  for  both  the  case  of  cre¬ 
ation  by  field  emission  and  of  creation  by  impact  ionization 
of  deep  levels.  Finally,  band-to-band  impact  ionization  is 
considered.  It  should  be  noted  that  both  interfaces  and  the 
space-charge  layer  are  assumed  to  be  discrete  traps. 

A.  Insulator-semiconductor  interfaces 

Emission  rates  from  the  insulator-phosphor  interface  are 
controlled  by  both  the  field  in  the  phosphor  layer  adjacent  to 
the  interface  of  interest  and  the  ambient  temperature.  At  low 
phosphor  fields,  phonon-assisted  tunneling  (PAT)  and  ther¬ 
mal  emission  (thermal)  of  carriers  may  be  important  mecha¬ 
nisms  for  carrier  emission.7  At  high  fields,  the  dominant 
mechanism  for  electron  emission  is  by  pure  tunnel  (PT)  in¬ 
jection  of  electrons  from  interface  states  into  the  conduction 
band  of  the  phosphor.  The  total  emission  rate  at  an  insulator- 
phosphor  interface  is  the  sum  of  the  three  separate  emission 
rates: 

(i§) 

Computation  of  the  separate  emission  rates  for  thermal  emis¬ 
sion,  phonon-assisted  tunneling,  and  pure  tunneling  is  per¬ 
formed  using  analytical  forms  reported  in  the  literature.20-22 
Then  with  the  total  emission  rate  established,  the  time  rate  of 
change  of  the  charge  on  an  interface  is  given  by 

-^f  =  -e„[qN0f0-qx(t)],  (19) 

where  NJ 0  is  the  no-field  occupancy  of  the  interface  in 
question,  qx(t)  is  the  charge  on  interface  x  at  time  t ,  and  en 
is  the  emission  rate  computed  from  Eq.  (18).  The  expressions 
employed  for  computation  of  ,  e^ermal,  and  e ^  are 
given  in  Appendix  A  along  with  a  derivation  of  Eq.  (19). 

B.  Space-charge  creation 

The  remaining  current  term  to  be  accounted  for  is  the 
space-charge  current  term,  dqjdt.  This  term  can  be  modeled 


in  one  of  two  ways  depending  on  whether  it  is  assumed  that 
space  charge  is  created  by  impact  ionization  of  deep-level 
traps  or  by  field  emission  from  states  in  the  band  gap  of  the 
phosphor  layer. 


1.  Space-charge  creation  by  impact  ionization 

Impact  ionization  processes  in  a  semiconductor  are  gen¬ 
erally  described  by  a  field-dependent  ionization  function  rep¬ 
resenting  the  number  of  impact  ionization  events  per  unit 
length.  To  achieve  a  more  practical  measure  of  impact  ion¬ 
ization  events  in  a  semiconductor  in  terms  of  externally  mea¬ 
sured  quantities,  the  ionization  function  is  manipulated  to 
obtain  a  multiplication  factor  across  the  region  of  interest. 
The  creation  of  space  charge  by  impact  ionization  of  deep- 
level  traps  is  described  by  an  ionization  function  of  the  form 


a, 


(/p)=kr 


trapW  p 


nT\  qfp 


N'l  Ei0 


exp 


n- 1,2, 


(20) 


where  nT  is  the  density  of  traps  containing  electrons,  Nt  is 
the  total  density  of  traps  in  the  phosphor,  £ion  is  an  effective 
ionization  energy,  and  f0  is  the  characteristic  field.  For  simu¬ 
lation  of  the  single-sheet  charge  model,  an  ionization  func¬ 
tion  is  computed  for  both  regions  of  the  phosphor,  and  the 
relevant  field,  fpl  or  fp2,  is  inserted  into  Eq.  (20).  Then, 
using  the  two  ionization  functions  computed  from  Eq.  (20),  a 
multiplication  factor  across  the  phosphor  layer  is  defined  as 


Mtrap(/pi  ,fp2)  =  exp[  ot^pifp  \)ds+  a^p(fp2)  (dp-ds)]. 

(21) 

Finally,  the  time  rate  of  change  of  space  charge,  dqjdt ,  is 
computed  from 


^=Jn(0,t)[Mtr3p(fpUfp2)-ll  (22) 

where  Jn(Q,t)  is  the  electron  current  emission  from  the  ca¬ 
thodic  interface.  It  should  be  noted  that  the  derivation  of  Eqs. 
(20)  through  (22)  is  presented  in  Appendix  B. 

In  summary,  Eqs.  (20)  through  (22)  are  used  to  model 
deep-level  trap  impact  ionization.  Note  that  the  amount  of 
deep-level  impact  ionization  is  determined  by  assuming  that 
it  occurs  throughout  the  phosphor.  However,  once  the 
amount  of  trap  ionization  is  determined,  it  is  assumed  to 
occur  exclusively  at  the  single  sheet  of  charge.  This  method¬ 
ology  is  consistent  with  the  philosophy  of  the  single-sheet 
charge  model,  which  implicitly  assumes  that  a  charge  distri¬ 
bution  may  be  approximated  as  a  single  sheet  which  repre¬ 
sents  the  charge  centroid  of  the  actual  charge  distribution. 


2.  Space-charge  creation  by  field  emission 

The  second  mechanism  for  space-charge  creation,  field 
emission  from  traps  in  the  bulk  phosphor,  is  described  in  an 
analogous  manner  to  interface  emission.  The  general  equa¬ 
tions  for  field  emission  from  a  space-charge  layer  are  identi¬ 
cal  to  those  for  field  emission  from  an  interface.  However, 
the  depth  of  the  bulk  traps  and  the  number  of  trapped  elec¬ 
trons  is  not,  in  general,  equal  to  the  values  for  the  interfaces. 
Therefore,  the  space-charge  layer  is  simply  treated  as  a  third 
interface  located  in  the  bulk  phosphor  layer. 

Keir,  Ang,  and  Wager  4671 


J.  Appl.  Phys.,  Vol.  78,  No.  7/1  October  1995 


C.  Band-to-band  impact  ionization 


In  the  range  of  phosphor  fields  at  which  ACTFEL  de¬ 
vices  operate,  electron  multiplication  by  band-to-band  im¬ 
pact  ionization  is  a  distinct  possibility.  An  often-used  form 
for  the  ionization  function,4 


/  r  \  qfp 

*band(//J  =  exp| 
^ion 


fo\ 
fp I 


n~  1,2, 


(23) 


is  employed  for  this  simulation.  In  the  above  equation,  Eion  is 
an  effective  ionization  energy  and  /0  is  a  constant  that  sets 
the  field  range  where  band-to-band  impact  ionization  be¬ 
comes  important.  Once  again,  since  the  focus  of  this  paper  is 
phosphor  space  charge,  the  multiplication  factor  associated 
with  band-to-band  impact  ionization  cannot  be  computed  us¬ 
ing  a  constant  field  assumption.  Therefore,  the  formulation  of 
Eq.  (21)  is  used  for  the  multiplication  factor  due  to  band-to- 
band  impact  ionization  with  the  terms  replaced  by  aband 
terms.  Finally,  the  additional  currents  at  the  interfaces  due  to 
band-to-band  impaction  ionization  are  computed  similar  to 
the  time  rate  of  change  of  space  charge  in  Eq.  (22),  except 
that  the  band-to-band  multiplication  factor  is  used.  It  should 
also  be  noted  that  the  holes  generated  by  the  band-to-band 
process  are  assumed  to  drift  to  the  cathodic  interface  and 
instantaneously  recombine  with  electrons  residing  at  the  ca¬ 
thodic  interface  in  this  simulation.  This  assumption  of  instan¬ 
taneous  electron-hole  recombination  is  not  appropriate  for 
simulation  of  maximum  charge-maximum  voltage 
(Gmax'^max)  characteristics4,23  but  is  appropriate  for  the  cur¬ 
rent  purpose  of  simulating  C-V  and  Q-Fp  overshoot  since 
recombination  does  not  significantly  affect  overshoot. 


IV.  FEEDBACK  MECHANISMS  OF  CARRIER  EMISSION 

The  carrier  emission  mechanisms  previously  discussed 
exhibit  a  feedback  effect  such  that  after  carriers  are  emitted, 
future  emission  rates  are  affected.  Each  of  the  methods  of 
carrier  emission,  interface  emission,  field  emission  from  bulk 
traps,  impact  ionization  of  bulk  traps,  and  band-to-band  im¬ 
pact  ionization  shows  feedback  effects  of  varying  degree. 

A.  Interface  emission 

It  has  been  established  previously  that  as  an  ACTFEL 
device  begins  conducting  charge,  the  charge  that  piles  up  on 
the  interfaces  sets  up  a  counterfield  that  opposes  the  applied 
field.3  This  effect  limits  the  field  in  the  phosphor  layer  of  the 
ACTFEL  device  and  can  give  rise  to  the  so-called  field¬ 
clamping  effect  seen  in  devices  not  exhibiting  space-charge 
effects.  A  pictorial  representation  of  the  feedback  effect  of 
interface  emission  is  illustrated  in  Fig.  5.  As  shown  in  Fig.  5, 
the  applied  voltage  vg  results  in  a  field  across  the  phosphor 
layer  which  causes  electron  emission  from  insulator- 
phosphor  interface  1.  Electrons  emitted  from  interface  1  drift 
across  the  phosphor  layer  and  are  collected  at  interface  2. 
Continued  emission  from  interface  1  leads  to  a  net  negative 
charge  on  interface  2,  and  a  net  positive  charge  on  interface 
1.  The  charge  present  at  the  two  interfaces  gives  rise  to  the 
counterfield  that  reduces  the  net  field  across  the  entire  phos¬ 
phor  layer.  This  self-limiting  effect  concomitant  with  inter¬ 


FIG.  5.  Interface  feedback  in  ACTFEL  devices. 


face  emission  can  be  viewed  as  a  negative  feedback  effect 
and  is  the  essence  of  the  feedback  mechanism  for  interface 
emission. 

B.  Space-charge  creation 

The  effects  of  the  presence  of  a  sheet  of  space  charge  on 
the  phosphor  field  can  be  seen  in  Fig.  6;  the  field  is  no  longer 
constant  across  the  phosphor  layer  as  it  is  in  the  absence  of 
space  charge;  rather,  the  phosphor  field  changes  discontinu- 
ously  across  the  sheet  of  space  charge.  This  discontinuous 
change  in  the  phosphor  field  across  the  sheet  of  charge  ap¬ 
plies  to  the  creation  of  space  charge  either  by  impact  ioniza¬ 
tion  of  deep-level  traps  or  by  field  emission  from  bulk  traps. 
Differences  between  the  two  mechanisms  arise  due  to  the 
manner  in  which  each  space-charge  creation  mechanism  re¬ 
acts  to  changes  in  phosphor  field  due  to  the  presence  of  bulk 
space  charge,  as  discussed  in  the  following  sections. 


FIG.  6.  Phosphor  field  perturbations  caused  by  the  presence  of  space  charge. 
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1 .  Impact  ionization  of  deep  levels 

The  creation  of  space  charge  by  impact  ionization  of 
deep-level  traps  has  the  weakest  feedback  effect  of  the  four 
mechanisms  of  carrier  emission  discussed.  The  feedback  ef¬ 
fects  of  carrier  emission  as  space  charge  is  created  by  impact 
ionization  of  deep  levels  is  best  described  by  first  consider¬ 
ing  the  effect  that  space  charge  has  on  the  phosphor  field 
components  and  then  proceeding  from  there.  As  space  charge 
is  created  by  impact  ionization  of  deep  levels,  the  field  be¬ 
tween  the  cathodic  interface  and  the  sheet  of  space  charge  is 
increased,  while  the  phosphor  field  between  the  sheet  of 
space  charge  and  the  anodic  interface  is  reduced.  These 
phosphor  field  perturbations  due  to  the  creation  of  space 
charge  lead  to  increased  electron  emission  from  and  trap  ion¬ 
ization  near  the  cathodic  interface.  Similarly,  the  reduced 
field  near  the  anodic  interface  results  in  decreased  deep-level 
impact  ionization  in  this  low-field  region.  Overall,  however, 
the  creation  of  space  charge  by  impact  ionization  of  deep 
levels  is  not  a  self-limiting  process  in  terms  of  field  feed¬ 
back.  In  fact,  the  feedback  in  this  situation  can  be  positive 
due  to  the  increased  interface  emission  and  the  increased 
ionization  rate  in  the  phosphor  region  near  the  cathodic  in¬ 
terface.  Fortunately,  the  interface  feedback  effect  helps  this 
process  from  running  away  by  reduction  of  the  field  across 
the  entire  phosphor  with  the  pileup  of  charge  at  the  inter¬ 
faces.  However,  the  effect  that  limits  the  creation  of  space 
charge  by  impact  ionization  of  deep  levels  is  the  finite  num¬ 
ber  of  traps  available  for  ionization  within  the  phosphor 
layer. 

2.  Field  emission  from  bulk  traps 

When  space  charge  is  created  by  field  emission  from 
bulk  traps,  the  phosphor  fields  on  either  side  of  the  space- 
charge  layer  are  affected  in  a  manner  similar  to  space-charge 
creation  by  impact  ionization  of  deep  levels.  The  similarity 
arises  because  the  creation  of  space  charge  increases  the 
electric  field  between  the  cathodic  interface  and  the  space- 
charge  layer  while  the  field  between  the  space-charge  layer 
and  the  anodic  interface  is  reduced.  When  space  charge  is 
created  by  field  emission  from  bulk  traps,  bulk  emission  of 
electrons  is  dependent  only  on  the  field  between  the  space- 
charge  layer  and  the  anodic  interface.  This  gives  rise  to  a 
negative  feedback  effect  when  space  charge  is  created  by 
field  emission  because  the  emission  of  electrons  from  the 
space-charge  layer  reduces  the  field  that  future  emission 
from  the  space-charge  layer  is  dependent  on.  Also,  the  in¬ 
creased  field  near  the  cathodic  interface  causes  increased 
emission  from  the  cathodic  interface,  which  consequently  in¬ 
creases  the  counterfield  across  the  entire  device,  as  discussed 
in  the  interface  emission  feedback  section. 

C.  Band-to-band  impact  ionization 

The  mechanism  with  the  strongest  feedback  effect  is 
band-to-band  impact  ionization.  The  exponential  rise  of  the 
multiplication  factor  as  a  function  of  field  can  provide  a  very 
large  electron  current  to  the  anodic  interface  along  with  a 
large  hole  current  to  the  cathodic  interface.  The  operation  of 
the  feedback  effect  associated  with  band-to-band  impact  ion¬ 
ization  is  similar  to  the  feedback  effect  of  interface  emission 


in  that  the  carriers  that  are  emitted  produce  a  counterfield 
across  the  entire  thickness  of  the  phosphor  layer  opposing 
the  applied  field.  However,  the  large  multiplication  associ¬ 
ated  with  band-to-band  impact  ionization  can  provide  many 
times  the  number  of  carriers  emitted  from  the  interfaces  and, 
hence,  a  much  stronger  feedback  effect  than  interface  emis¬ 
sion  due  to  the  sheer  number  of  carriers  involved.  Further¬ 
more,  the  band-to-band  impact  ionization  mechanism  can  be 
a  major  factor  in  keeping  space-charge  creation  by  impact 
ionization  of  deep  levels  at  a  steady  level. 


V.  ELECTRICAL  CHARACTERIZATION  OF  ACTFEL 
DEVICES 


The  focus  of  this  paper  is  the  inclusion  of  bulk  space- 
charge  generation  into  the  standard  ACTFEL  device  model 
to  assess  its  impact  on  a  device’s  C-V  and  Q-Fp  character¬ 
istics.  However,  several  issues  need  to  be  clarified  about  the 
actual  C-V  and  Q-Fp  measurements  themselves  in  relation 
to  space  charge. 

A.  C-  V  analysis  with  space  charge 

A  C-V  measurement  is  performed  by  monitoring  the 
current  through  the  ACTFEL  device  during  the  rising  edge  of 
the  applied  voltage  waveform.  Then,  the  measured  capaci¬ 
tance  of  the  ACTFEL  device  is  determined14  by  the  relation 


*(*) 

dvgldt * 


(24) 


However,  the  quantity  i(t)  is  known  in  terms  of  the  single¬ 
sheet  charge  model  as  dqjdt.  Therefore,  when  Eq.  (12)  is 
inserted  into  Eq.  (24),  an  expression  for  the  measured  device 
capacitance  in  terms  of  the  single-sheet  charge  model 
emerges, 


r  — 


ciF  cp 


dq\ 

dVo 


+  1- 


s< 

dvQ 


+  Ct. 


(25) 


The  above  equation  shows  that  space-charge  creation  adds  to 
the  measured  capacitance  through  the  dqjd vg  term  and  is 
the  mechanism  responsible  for  C-V  overshoot.  Note  that  Eq. 
(25)  indicates  that  the  location  where  the  space  charge  is 
generated  also  determines  the  magnitude  of  the  C-V  over¬ 
shoot;  the  maximum  C-V  overshoot  will  occur  when  the 
sheet  of  charge  is  located  near  the  cathodic  interface. 


B.  Q-Fp  analysis  with  space  charge 

Recently,  the  Q-Fp  technique  has  been  demonstrated  to 
be  a  useful  electrical  measurement  for  ACTFEL  device 
characterization.13-15,24  However,  Q-Fp  analysis  implicitly 
assumes  that  no  space  charge  exists  in  the  ACTFEL  phos¬ 
phor  layer  and,  thus,  that  the  electric  field  is  constant  across 
the  phosphor  layer.3  Since  this  is  not  always  a  reasonable 
assumption,  it  needs  to  be  determined  what  is  being  mea¬ 
sured  by  this  technique  when  space  charge  is  present  in  the 
phosphor.  Thus  the  purpose  of  this  section  of  the  paper  is  to 
determine  what  qmt(t)  and  fp(t)  in  a  Q-Fp  measurement 
actually  represent  when  space  charge  is  present  in  the  phos¬ 
phor  layer. 
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A  Q-Fp  curve  is  created  by  plotting  the  following  ex¬ 
pression: 


C/  +  C0 


qe(t)-c  vg(t) 


From  Eqs.  (4)  and  (5),  it  is  determined  that 

cpdpfpi{t)=-q2{t)-qe(t).  (34) 

Then,  substitution  of  the  above  equation  and  Eq.  (3)  into  Eq. 
(32)  yields 


versus 


1  qe(t) 

fp(t)=d~  ~c - Vs(t)  ■  (27) 

up  l  ci 

From  the  above  equations,  it  is  seen  that  the  measurable 
quantities  qe(t )  and  vg(t)  are  needed  to  compute  the  quan¬ 
tities  q int( t )  Therefore,  to  determine  what  is  being 

measured  when  these  quantities  are  plotted  against  one  an¬ 
other  (within  the  framework  of  the  single-sheet  charge 
model),  it  is  necessary  to  solve  Eqs.  (1)  through  (6)  for  the 
quantities  on  the  right-hand  sides  of  Eqs.  (26)  and  (27). 

To  accomplish  this,  the  right-hand  side  of  Eq.  (27)  is  first 
solved  to  determine  what  fp(t)  represents  in  a  device  con¬ 
taining  a  sheet  of  space  charge.  For  an  ACTFEL  device  with 
identical  insulating  layers,  it  is  known  that 

cndnfn(t)  =  cadi2fi2(t)  =  -qe(t).  (28) 

For  an  ACTFEL  device,  the  series  combination  of  the  insu¬ 
lator  capacitances  leads  to  the  total  insulator  capacitance. 
Therefore,  for  an  ACTFEL  device  with  identical  insulating 
layers, 

Cii  —  Cft—lCi ,  (29) 

where  ct  is  the  insulator  capacitance.  Plugging  the  above 

equation  into  Eq.  (7),  it  is  found  that 

dnfn(t)=di2fi2(t)  =  -^-.  (30) 

Then,  the  above  equation  is  substituted  into  Eq.  (1)  to  yield 

-^~+dJpl(t)+fp2(t)(dp-ds)=-vg(t).  (31) 

c  i 

Rearranging  the  above  equation,  it  is  found  that 


dP 


fp\(*)  + 


1 

dp 


(32) 


as  desired.  Comparing  the  right-hand  side  of  Eq.  (32)  to  Eq. 
(27),  it  is  seen  that  the  actual  quantity  measured  for  the  phos¬ 
phor  field  when  a  Q-Fp  measurement  is  performed  is  the 
average  phosphor  field.  This  result  is  derived  from  the 
single-sheet  charge  model,  but  it  can  be  shown  inductively 
that  for  n  layers  of  sheet  charge  the  result  obtained  in  the 
average  field  is  still  what  is  measured. 

Next,  it  is  necessary  to  determine  what  the  qint  term  in 
the  Q-Fp  plot  represents  when  space  charge  is  present  in  the 
phosphor  layer.  To  determine  this  in  terms  of  the  single-sheet 
charge  model,  it  is  necessary  to  solve  for  the  right-hand  side 
of  Eq.  (26)  using  Eqs.  (1)  through  (6).  Rearranging  Eq.  (32), 
it  is  found  that 


^  pd  pf pi(t) 


Cp<le(t) 


cpVg(t)  cpds\.f p\(t)  fp2(t)^* 

(33) 


<lsc(‘)=  ~  qe(t)-cpvg(t).  (35) 

ap  W 

Additionally,  substitution  of  the  integral  of  Eq.  (12)  into  the 
above  equation  gives  qint(t)  in  terms  of  qx{t)  and  qsc(t)y  as 
follows: 

gl(0  +  |l  +  ^jgsc(0=  Pc  '  qe(t)-CpV  git)-  (36) 

Therefore,  a  comparison  of  Eqs.  (36)  and  (26)  shows  that  the 
internal  charge  measured  in  a  Q-Fp  plot  of  a  sample  exhib¬ 
iting  space-charge  effects  is  equal  to  the  actual  internal 
charge  plus  a  factor  which  depends  on  both  the  amount  and 
the  location  of  the  space  charge  in  the  phosphor  layer. 

VI.  SIMULATION  IMPLEMENTATION 

The  ACTFEL  device  modeling  simulation  results  are  ob¬ 
tained  from  a  realization  of  the  single-sheet  charge  model 
written  in  the  C  programming  language.  The  simulation  re¬ 
sults  are  obtained  from  a  numerical  solution  of  Eqs.  (13)  and 
(14)  with  a  fourth/fifth  order  Runge-Kutta-Fehlberg  algo¬ 
rithm.  All  simulations  are  written  and  performed  using  a 
Hewlett-Packard  9000  type  computer.  Simulation  of  ten  pe¬ 
riods  of  a  1  kHz  driving  waveform  takes  between  15  s  and  2 
min,  depending  on  the  program  being  run  and  system  usage. 

The  parameters  used  in  the  simulations  are  listed  in 
Table  I.  It  should  be  noted  that  the  parameters  Eion ,  E[on ,  Nt , 
F0,  and  F*0  apply  only  to  the  simulation  that  assumes  space- 
charge  creation  by  impact  ionization.  Similarly,  the  param¬ 
eters  NofQ(sc)  and  Esc  apply  only  to  the  simulation  that 
assumes  space  charge  creation  by  field  emission  from  bulk 
traps.  The  remaining  parameters  apply  to  both  types  of  simu¬ 
lation.  Finally,  the  ‘V’  term  in  Eqs.  (19)  and  (22)  is  assumed 
to  be  unity  for  fields  below  the  field  threshold,  and  two  for 
fields  above  this  threshold,25  where /0  in  Eqs.  (19)  and  (22) 
is  taken  to  be  the  threshold  field. 

The  simulation  is  performed  for  the  case  of  aluminum- 
titanium  oxide  (ATO)  insulating  layers,  and  a  ZnS:Mn  phos¬ 
phor  layer.  The  dielectric  constants  of  these  materials  are 
used  in  conjunction  with  the  thicknesses  of  each  layer  of  the 
ACTFEL  device  to  determine  the  capacitance  of  each  layer. 
The  thicknesses  of  each  layer  are  chosen  to  correspond  to  the 
devices  tested  by  Abu-Dayah  et  al 14 

VII.  RESULTS  AND  DISCUSSION 

A.  The  origins  of  C-  V  and  Q-Fp  overshoot 

One  of  the  primary  motivating  issues  for  undertaking  the 
work  presented  herein  is  the  determination  of  the  origins  of 
C-V  and  Q-Fp  overshoot.  First,  consider  aC-V  measure¬ 
ment  of  an  ACTFEL  device.  The  measurement  is  performed 
by  measuring  the  instantaneous  current  through  the  ACTFEL 
device  and  dividing  this  value  by  the  time  rate  of  change  of 
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TABLE  I.  Nominal  parameter  values  used  for  implementing  the  single-sheet  charge  model. 


Parameter 

Description 

Nominal  parameter  value 

da 

Thickness  of  insulator  1 

2700  A 

da 

Thickness  of  insulator  2 

2700  A 

dP 

Thickness  of  phosphor  layer 

5500  A 

ds 

Space-charge  location 

3500  A 

Nofoi'P) 

No-field  interface  occupation 

5X1013  cm-2 

N, 

Bulk  trap  concentration 

lX1017cnT3 

sccf 

Space-charge  capture  probability 

0.1 

£« 

Interface  trap  depth 

1.0  eV 

Wo/oCsc) 

No-field  space  charge  occupation 

1X1012  cm"2 

Space-charge  trap  depth 

0.9  eV 

^ion 

Band-to-band  effective  ionization  energy 

5.4  eV 

EL 

Trap-to-band  effective  ionization  energy 

2.6  eV 

F0 

Band-to-band  characteristic  field 

1.7  MV/cm 

F'o 

Trap-to-band  characteristic  field 

3.0  MV/cm 

the  voltage  dropped  across  the  ACTFEL  device.  When  space 
charge  is  being  created  in  the  ACTFEL  device  it  perturbs  the 
measured  current  as  seen  in  Eqs.  (12)  and  (25),  so  that  a 
capacitance  greater  than  the  insulator  capacitance  may  be 
observed  near  turn-on.  However,  as  discussed  in  the  section 
on  feedback  effects  of  space  charge,  space-charge  emission 
occurs  such  that  the  creation  of  space  charge  tends  to  turn 
itself  off.  Therefore,  as  the  rate  of  space-charge  creation  falls 
towards  zero  above  turn-on,  the  measured  current  falls  to  a 
level  close  to  what  would  be  measured  if  no  space  charge  is 
present.  Hence,  this  produces  overshoot  in  the  C-  V  charac¬ 
teristic  of  an  ACTFEL  device  with  space  charge  present  in 
the  phosphor. 

The  issue  of  field  overshoot  in  Q-Fp  plots  presents  an¬ 
other  puzzling  problem  without  the  inclusion  of  space-charge 
effects.  The  origins  of  field  overshoot  in  Q-Fp  plots  is  simi¬ 
lar  in  nature  to  the  overshoot  in  C-V  plots,  but  relates  to  a 
different  measured  quantity.  The  field  overshoot  in  a  Q-Fp 
plot  arises  due  to  the  fact  that  the  Fp  measured  in  a  Q-Fp 
plot  is  actually  the  average  phosphor  field.  The  reason  that 
overshoot  is  seen  is  that  the  average  phosphor  field  is  at  first 
a  certain  value  that  is  then  reduced  after  space-charge  cre¬ 
ation  begins.  This  can  again  be  traced  to  the  feedback  effect 
of  space-charge  creation. 

B.  Impact  ionization 

Simulation  of  space  charge  in  the  phosphor  layer  by  im¬ 
pact  ionization  of  deep-level  impurities  shows  several  inter¬ 
esting  effects.  First,  the  addition  of  band-to-band  impact  ion¬ 
ization  to  the  simulation  along  with  impact  ionization  of 
traps  leads  to  a  significantly  larger  overshoot  in  both  the  C-V 
and  Q-Fp  plots  than  when  only  impact  ionization  of  traps  is 
included.  A  comparison  of  simulated  C-V  plots  both  with 
and  without  band-to-band  impact  ionization  is  given  in  Fig. 
7. 

It  is  established  in  Eq.  (24)  that  the  magnitude  of  the 
C-V  overshoot  is  a  function  of  the  location  of  the  space- 
charge  layer  and  the  partial  derivative  of  the  amount  of  space 
charge  with  respect  to  the  applied  voltage.  In  steady  state,  the 
amount  of  space  charge  generated  during  one  period  of  the 
driving  waveform  is  equal  to  the  amount  of  space  charge 
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annihilated  during  the  period.  Therefore,  when  band-to-band 
impact  ionization  is  included  in  the  simulation  along  with 
impact  ionization  of  traps,  there  is  a  significantly  larger 
amount  of  charge  transferred  to  the  insulator-phosphor  inter¬ 
faces.  This  increase  in  charge  on  the  interfaces  leads  to  a 
larger  amount  of  charge  transferred  during  the  interpulse  in¬ 
terval  or  low-field  portion  of  the  driving  waveform.  The 
larger  amount  of  charge  transfer  during  the  interpulse  inter¬ 
val  leads  to  a  greater  amount  of  space-charge  annihilation, 
and  hence,  a  greater  overshoot  because  the  amount  of  space 
charge  annihilated  in  one  period  must  equal  the  amount  of 
space  charge  generated  in  steady  state. 

A  simulated  Q-Fp  plot  with  both  band-to-band  and 
deep-level  impact  ionization  included  and  a  simulated  Q-Fp 
plot  without  bulk  space  charge  are  shown  in  Fig.  8.  Figure  8 
shows  that  space-charge  creation  by  impact  ionization  causes 
the  average  phosphor  field  for  different  polarity  pulses  to  be 
quite  different.  In  contrast,  the  simulated  Q-Fp  curve  with¬ 
out  space  charge  included  in  the  phosphor  layer  shows  that 
the  field  for  both  polarity  pulses  is  approximately  equal. 
Also,  Fig.  8  shows  that  more  charge  is  transferred  across  the 
phosphor  layer  when  space  charge  is  present,  as  would  be 
expected  intuitively.  Finally,  Fig.  8  shows  that  the  Q-Fp 
curve  is  shifted  so  it  is  no  longer  centered  around  the  point 


FIG.  7.  Simulated  Al-  C-V  plots  of  an  ACTFEL  device  with  (a)  both 
impact  ionization  of  deep-level  traps  and  band-to-band  impact  ionization 
included  and  (b)  only  impact  ionization  of  deep-level  traps  included. 
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FIG.  8.  Simulated  Q-F p  plots  of  an  ACTFEL  device  with  (a)  both  impact 
ionization  of  deep-level  traps  and  band-to-band  impact  ionization  included 
and  (b)  impact  ionization  of  traps  not  included. 

where  the  internal  charge  and  the  phosphor  field  are  zero 
when  space-charge  creation  is  included. 

C.  Field  emission 

The  creation  of  space  charge  by  field  emission  from  bulk 
traps  is  the  second  type  of  single-sheet  charge  model  simu¬ 
lated.  As  is  the  case  of  deep  level  plus  band-to-band  impact 
ionization,  the  creation  of  space  charge  by  field  emission 
causes  both  C-V  and  Q-F p  overshoot.  However,  for  space- 
charge  creation  by  field  emission  it  is  not  necessary  to  have 
band-to-band  impact  ionization  present  in  order  to  observe 
field  overshoot.  This  is  due  to  the  fact  that  space-charge  cre¬ 
ation  by  field  emission  effectively  limits  itself,  as  discussed 
in  the  section  on  feedback  effects  of  space-charge  creation. 

The  simulated  results  for  space-charge  creation  by  field 
emission  are  very  similar  to  those  simulated  assuming  im¬ 
pact  ionization;  this  can  be  seen  by  comparing  the  C-V  plots 
of  space-charge  creation  by  field  emission  shown  in  Fig.  9  to 
the  corresponding  plots  for  impact  ionization  in  Fig.  7.  It 
should  be  noted  that  the  device  structures  for  both  the  field 
emission  and  the  impact  ionization  simulations  shown  in 
Figs.  7  and  9  are  identical.  Furthermore,  the  space-charge 
layer  is  located  at  the  same  position  in  the  phosphor  layer  for 
both  simulations.  For  comparison  purposes,  a  simulated  C-V 
plot  without  any  space  charge  is  shown  in  Fig.  10.  The  rel¬ 
evance  of  this  plot  is  that  exactly  the  same  parameters  are 


FIG.  9.  Simulated  C-  V  plots  of  both  polarity  pulses  for  an  ACTFEL  device 
with  a  space-charge  layer  at  ds= 3500  A  which  is  created  by  field  emission. 


FIG.  10.  Simulated  C-  V  plot  of  an  ACTFEL  device  without  space  charge  in 
the  phosphor  layer. 

used  for  the  simulations  shown  in  Figs.  9  and  10.  However, 
as  seen,  the  presence  of  space  charge  in  the  phosphor  leads  to 
a  C-V  response  exhibiting  overshoot. 

A  Q-F p  plot  with  bulk  space  charge  created  by  field 
emission  is  shown  in  Fig.  11.  It  is  seen  in  Fig.  11  that  there 
is  again  a  small  field  overshoot  for  one  field  polarity,  as  is  the 
case  for  space  charge  created  by  impact  ionization.  Further¬ 
more,  the  maximum  average  phosphor  fields  shown  in  the 
Q-F p  plot  are  again  different  for  each  polarity  of  field.  Also, 
it  is  seen  that  space-charge  creation  by  field  emission  makes 
the  Q-F p  plot  asymmetrical  about  the  line  where  internal 
charge  is  zero.  Furthermore,  it  is  seen  that  the  Q-F p  plot  of 
an  ACTFEL  device  without  space  charge  present  in  the  phos¬ 
phor  layer  exhibits  symmetry  about  the  line  where  the  inter¬ 
nal  charge  is  zero.  The  result  of  these  asymmetries  is  that  the 
Q-F p  plot  is  offset  so  that  the  center  of  the  Q-F p  curve  is  no 
longer  the  point  where  both  the  internal  charge  and  the  phos¬ 
phor  field  are  zero,  as  it  is  in  the  Q-F p  curve  for  an  ACTFEL 
device  with  no  space  charge  in  the  phosphor.  It  is  interesting 
to  note  that  a  similar  offset  of  internal  charge  is  not  seen  in 
the  case  of  space-charge  creation  by  impact  ionization  of 
deep-level  traps,  as  evident  in  Fig.  8. 

D.  Space-charge  parameter  variation  trends 

As  mentioned  previously,  the  parameters  used  in  the 
simulation  of  Figs.  7  and  10  are  collected  in  Table  I.  The 
purpose  of  the  following  subsection  is  to  explore  simulation 


FIG.  11.  Simulated  Q-F p  plots  of  an  ACTFEL  device  with  (a)  a  layer  of 
space  charge  at  c^=3500  A  created  by  field  emission  and  (b)  no  space 
charge  in  the  phosphor  layer. 
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FIG.  12.  Family  of  C-V  curves  for  Al-  (the  arrow  indicates  increasing  ds ; 
^=500,  1000,  1500,  2750,  4000,  4500,  and  5000  A). 

trends  when  the  magnitudes  of  some  of  the  key  simulation 
parameters  are  systematically  varied.  It  should  be  noted  that 
the  following  discussions  of  variations  of  space-charge  pa¬ 
rameters  assume  the  device  structure  of  Fig.  1,  where  Vg  is 
applied  to  an  A1  electrode  and  the  indium  tin  oxide  (ITO) 
electrode  is  grounded. 

1 .  Trends  when  ds  is  varied 

The  results  of  a  parametric  variation  of  ds ,  the  location 
of  the  space-charge  layer,  on  the  C-  V  plots  of  both  polarities 
is  provided  in  Figs.  12  and  13.  As  seen  in  these  figures,  the 
location  of  the  space  charge  has  a  very  important  effect  on 
the  magnitude  of  the  overshoot  seen  in  the  C-V  plot.  Also, 
the  location  of  the  space  charge  determines  the  polarity  of 
the  pulse  for  which  overshoot  is  witnessed.  For  example,  the 
family  of  C-V  plots  shown  in  Fig.  12  demonstrates  that  the 
overshoot  of  a  negative  polarity  pulse  increases  as  the  space- 
charge  layer  nears  the  insulator-phosphor  interface  on  the 
ITO  side  of  the  device.  Similarly,  Fig.  13  shows  that  the 
overshoot  exhibited  in  a  C-V  curve  from  a  positive  applied 
pulse  decreases  as  the  space-charge  centroid  moves  closer  to 
the  insulator-  phosphor  interface  on  the  ITO  side  of  the  de¬ 
vice.  This  trend  is  actually  predicted  by  Eq.  (25)  which 
shows  that  the  C-  V  overshoot  term  is  actually  the  product  of 
dqJdVg  and  the  location  dependent  term,  1 —(dsldp). 

A  similar  trend  is  apparent  for  the  family  of  Q-Fp  plots 
shown  in  Fig.  14.  As  the  space-charge  layer  moves  close  to 
one  interface,  the  overshoot  observed  is  very  asymmetrical  in 


FIG.  13.  Family  of  C-V  curves  for  A1+  (the  arrow  indicates  increasing  ds ; 
dj=500,  1000,  1500,  2750,  4000,  4500,  and  5000  A). 
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FIG.  14.  Family  of  Q-Fp  curves  representing  changing  ds  (the  arrow  indi¬ 
cates  increasing  ds ;  ^=500,  2750,  and  5000  A). 

that  it  is  seen  for  only  one  polarity  pulse.  This  asymmetry 
appears  near  both  the  far  left  and  right  side  of  the  Q-Fp  plot 
and  the  arrows  mark  increasing  ds .  Therefore,  for  a  labora¬ 
tory  measurement,  the  location  of  the  space  charge  in  the 
phosphor  layer  can  be  roughly  determined  by  examining  the 
asymmetry  in  the  overshoot  characteristics  and  comparing  it 
to  simulated  results. 

Comparison  of  the  simulated  results  shown  in  Figs.  12 
and  13  with  the  results  from  experiments  performed13,14  on 
ALE  ZnS:Mn  devices  by  Abu-Dayah  et  al.  shed  some  light 
on  the  location  of  the  space-charge  layer  in  these  devices. 
Figure  13  shows  that  small  or  no  overshoot  is  seen  for  an 
A1+  C-V  measurement  when  the  space-charge  layer  is  lo¬ 
cated  near  the  insulator-phosphor  interface  on  the  Al  side  of 
the  device.  Alternately,  Fig.  12  shows  that  a  large  overshoot 
is  produced  for  an  Al—  C-V  measurement  when  the  space- 
charge  layer  is  located  near  the  insulator-phosphor  interface 
on  the  ITO  side  of  the  ACTFEL  device.  When  Figs.  12  and 
13  are  compared  with  the  experimental  results  of  Abu-Dayah 
et  a/.,13  it  is  seen  that  the  experimental  results  are  very  simi¬ 
lar  to  the  simulated  results  when  the  sheet  of  space  charge  is 
located  very  near  the  insulator-phosphor  interface  on  the  Al 
side  of  the  device.  Thus,  the  simulated  results  provide  evi¬ 
dence  that  the  centroid  of  the  space-charge  distribution  for 
the  ACTFEL  device  samples  used  for  the  study  by  Abu- 
Dayah  et  al.  lies  near  the  insulator-phosphor  interface  on  the 
Al  side  of  the  device;  this  result  is  in  contrast  to  previous 
conclusions  of  Douglas  et  al.26  and  Abu-Dayah  et  alnM 
who  asserted  that  space-charge  creation  occurs  near  the  bot¬ 
tom,  ITO  interface. 

2.  Trends  exhibited  with  variation  of  capture 
efficiency 

Another  space-charge  simulation  parameter  is  the  space- 
charge  capture  efficiency.  Presently,  the  simulation  assumes 
that  a  negligible  amount  of  space  charge  is  annihilated  dur¬ 
ing  the  application  of  voltage  pulses  to  the  ACTFEL  device. 
However,  in  the  interpulse  interval  it  is  assumed  that  space 
charge  is  likely  to  be  annihilated  and  can  be  described  by  a 
capture  ratio  that  determines  what  percentage  of  the  charge 
transferred  during  the  interpulse  interval  annihilates  space 
charge.  This  assumption  is  made  because  prior  researchers 
have  shown  that  the  capture  cross  section  for  space-charge 
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FIG.  15.  Family  of  Q-Fp  curves  representing  a  variation  of  the  bulk  trap 
capture  efficiency  (the  arrows  indicate  increasing  capture  efficiency;  sccf 
=0.01,0.05,0.1,0.3,0.5,0.7). 


capture  processes  decreases  rapidly  with  field.27  Therefore, 
during  the  high-field  portion  of  the  operation  of  the  ACTFEL 
device  it  is  assumed  that  the  amount  of  electron  capture  by 
the  space-charge  layer  is  negligible. 

According  to  simulated  results  this  capture  ratio  for 
space-charge  annihilation  is  a  very  important  parameter.  The 
results  of  several  simulations  which  show  the  effects  of  this 
parameter  are  illustrated  in  Fig.  15.  These  simulations  show 
that  as  space-charge  capture  becomes  more  efficient,  the  de¬ 
vice  exhibits  increasing  overshoot.  Also,  as  capture  becomes 
more  efficient,  the  interpulse  leakage  charge  is  reduced.  The 
basic  reasoning  behind  this  effect  is  once  again  related  to  the 
amount  of  space  charge  being  constant  at  points  in  time  sepa¬ 
rated  by  one  period  of  the  driving  waveform.  This  means  that 
as  the  space  charge  becomes  more  efficient  at  charge  capture 
during  the  low-field  portion  of  the  pulse,  more  space  charge 
that  will  be  annihilated.  However,  because  the  device  is  in 
steady  state,  the  amount  of  space  charge  annihilated  during 
the  interpulse  interval  must  equal  the  amount  of  space  charge 
created  during  the  pulse.  Therefore,  greater  overshoot  is  wit¬ 
nessed  in  accordance  with  Eq.  (25).  This  result  leads  to  the 
conclusion  that  it  is  not  the  presence  of  space  charge  in  the 
phosphor  layer  that  causes  many  of  the  effects  attributed  to 
space  charge,  but  the  interaction  between  the  charge  at  the 
sheet  of  space  charge  and  at  the  insulator-phosphor  inter¬ 
faces. 


VIII.  CONCLUSIONS 

Simulation  of  space  charge  in  the  phosphor  layer  of  an 
ACTFEL  device  using  a  single-sheet  charge  model  provides 
evidence  that  space-charge  creation  is  indeed  responsible  for 
C-V  and  Q-Fp  overshoot.  Specifically,  Figs.  7  and  10  show 
that  space-charge  creation  by  either  impact  ionization  of 
deep-level  traps  or  by  field  emission  from  bulk  traps  causes 
overshoot  in  C-V  measurements  of  ACTFEL  devices.  In  ad¬ 
dition,  simulated  Q-Fp  overshoot  is  seen  in  Figs.  8,  11,  14, 
and  15.  Furthermore,  the  amount  of  overshoot  exhibited  by  a 
sample  is  dependent  on  several  parameters  such  as  the  loca¬ 
tion  of  the  space  charge  in  the  phosphor  layer  and  the  capture 
efficiency  of  the  bulk  traps.  The  amount  of  overshoot  exhib¬ 
ited  by  a  sample  also  depends  on  the  polarity  of  the  applied 


voltage  pulse.  As  shown  in  Fig.  10,  C-V  plots  for  different 
polarities  of  the  applied  voltage  can  exhibit  vastly  different 
overshoot  characteristics. 

The  simulation  of  space  charge  in  the  phosphor  layer 
using  a  single-sheet  charge  model  has  indicated  that  a  space- 
charge  distribution  with  a  centroid  near  one  of  the  insulator- 
phosphor  interfaces  causes  the  asymmetrical  overshoot  char¬ 
acteristics  witnessed  in  Fig.  9  and  experimental  work.13-15 
From  the  overshoot  characteristics  of  two  different  polarity 
C-  V  plots  for  a  given  ACTFEL  device  sample,  a  rough  esti¬ 
mate  of  the  location  of  the  centroid  of  the  space-charge  dis¬ 
tribution  in  the  phosphor  layer  can  be  determined  using  the 
simulated  data  given  in  Figs.  12  and  13.  Also,  the  amount  of 
overshoot  and  the  leakage  charge  exhibited  by  a  sample  has 
been  shown  to  be  dependent  on  the  capture  efficiency  of  the 
charged  bulk  phosphor  layer  traps.  However,  the  capture  ef¬ 
ficiency  does  not  show  any  effects  related  to  applied  voltage 
polarity,  as  does  the  location  of  the  space-charge  layer. 
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APPENDIX  A:  INTERFACE  EMISSION  RATES 


The  expression  for  the  emission  rate  due  to  pure  tunnel¬ 
ing  from  a  discrete  Coulombic  well  eY  is20“22 


PT/  r  \ 


4  (2 m*)mE%2 
3  qhfp 


X 


1- 


(A£it 

l  fit 


5/3 


(Al) 


where  Eit  is  the  depth  of  the  discrete  interfacial  trap  and 
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Next,  the  thermal  emission  rate  from  a  discrete  trap, 
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where  k  is  Boltzmann’s  constant,  T  is  the  temperature  in 
degrees  Kelvin,  a  is  the  capture  cross  section  for  thermal 
emission,  Nc  is  the  effective  density  of  states,  and  is  the 
thermal  velocity.  Finally,  the  emission  rate  due  to  phonon- 
assisted  tunneling  from  a  discrete  trap,  >  is  given  by20-22 
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These  emission  rates  given  by  Eqs.  (Al),  (A3),  and  (A4)  are 
then  inserted  into  Eq,  (18)  to  yield  the  total  emission  rate  of 
the  interface. 

The  time  rate  of  change  of  qx(t ),  the  charge  on  interface 
x  at  time  t ,  is 


dqx 

dt 


-qennx(t), 


(A5) 


where  nx{t)  is  the  number  of  electrons  present  at  interface  x 
at  time  t.  The  expression  for  the  charge  at  an  interface  in 
terms  of  n(t)  is  given  by3 

«I(0  =  «Wo-'>i(0].  (A6) 

where  Nq/q  is  the  no-field  occupancy  of  the  interface  in 
question.  Hence,  solving  for  nx(t)  in  the  above  equation  and 
substituting  into  Eq.  (A5)  yields 

^f  =  -ea[qN0f0-qx(t)],  (A7) 


where  qx(t)  can  be  found  in  terms  of  the  phosphor  fields 
using  Eqs.  (1)  through  (6). 


APPENDIX  B:  IMPACT  IONIZATION  OF  DEEP-LEVEL 
TRAPS 

The  impact  ionization  of  deep-level  traps  in  a  semicon¬ 
ductor  has  been  shown  to  be  fit  by  ionization  functions  of  the 
form25 

«trap(/p)  =  fi(/p)exp  — jjr)  .  n=  1,2.  (Bl) 

Also,  the  ionization  function,  atmp(fp)  is  often  expressed25,28 

<*traP  (fP)^Nta{fp)9  (B2) 

where  cr(fp)  is  the  capture  cross  section  for  the  process  and 
Nt  is  the  trap  density  in  the  semiconductor.  A  common  field- 
dependent  form  for  the  ionization  function  multiplier,  B(f  ), 
is4 

B{fp)= fk  (B3) 

^ion 

If  Eqs.  (Bl)  and  (B3)  are  used  to  determine  the  ionization 
function,  the  cross  section  for  the  process,  o{fp)  can  be 
found  by  equating  the  ionization  function  determined  by  Eqs. 
(Bl)  and  (B3)  to  that  in  Eq.  (B2).  This  results  in 

n=  1,2. 

(B4) 

However,  Eq.  (B4)  refers  to  a  situation  in  which  all  of  the 
traps  are  filled  with  electrons  and,  hence,  Nt  trapped  elec¬ 
trons  are  always  available  to  be  promoted  to  the  conduction 
band.  Since  it  is  assumed  that  the  impact  ionization  process 
creates  space  charge  in  the  phosphor  layer,  there  is  a  reduc¬ 
tion  in  the  number  of  traps  available  for  ionization  as  traps 
become  ionized.  Therefore,  the  quantity  nr,  the  density  of 
traps  containing  electrons,  must  be  introduced  into  Eq.  (B4) 
to  obtain  the  final  expression  for  the  impact  ionization  func¬ 
tion, 


(t  s  ...  /»r\  qfp  /of 
a^p(fP)  =  nTar(Jp)  =  \—\  —exp  -\—\  ,  n=  1,2. 

\lvtJ  ion  L  \Jpl  . 

(B5) 

With  the  form  of  the  ionization  function  now  estab¬ 
lished,  it  is  necessary  to  examine  current  continuity.  Neglect¬ 
ing  recombination-generation  processes  in  the  phosphor,  the 
current  continuity  equation  for  electrons  becomes4 

dJn 

=  Straps  0*6) 

which  has  the  solution 

Jn{x,t)  =  7„(0,r)exp[  (B7) 

However,  in  the  single-sheet  charge  model  it  is  implicitly 
assumed  that  there  are  two  regions  of  the  phosphor  layer 
whose  fields  are  not  in  general  equal.  Therefore,  it  is  neces¬ 
sary  to  examine  the  current  continuity  equation  in  each  re¬ 
gion  separately.  First,  in  region  p  1 , 

dJ  n 

•^-  =  atraps(/pl  )Jn,  0=Sx=£tf ,  (B8) 

meaning  that 

Jn{x,t)=Jn(0,t)exp[amps(fpi)x],  0^x*zds.  (B9) 
Similarly,  in  region  p2 , 
dJn 

=  «traps(/p2)^«  -  ds^X^dp  ,  (BIO) 

such  that 

Jn(x,t)  =  Jn(ds,t)exp[attaps{fp2)(x-ds)],  ds^x^dp. 

(Bll) 

Then,  since  Eqs.  (B9)  and  (Bll)  must  be  equal  at  ds,  it  is 
determined  that 

Jn(ds  ,t)  =  Jn{0,t)exp[aUzps(fpl)ds\.  (B12) 

Finally,  substituting  Eq.  (B12)  into  Eq.  (Bll)  and  evaluating 
the  resultant  expression  at  the  anodic  interface,  the  current 
arriving  at  the  anodic  interface  is  determined  to  be 

d n(dp  ,t)  —  ./fl(0,j)exp[ pi)d^\ 

x  exp[  atraps(/p2 )  (dp- ds)\.  (B13) 

Therefore,  the  multiplication  factor  across  the  entire  phos¬ 
phor  is  given  by 

p\  5/^2)  exp[<^traps(y*/?l)^j  atxapsif  pi) (dp  ds)~\ . 

(B14) 

It  should  be  noted  that  this  analysis  is  only  valid  for  the  case 
where  interface  1  is  the  cathodic  interface  and  interface  2  is 
the  anodic  interface  (i.e.,  a  positive  applied  voltage  to  the  Al 
contact).  However,  the  analysis  for  a  negative  applied  volt¬ 
age  proceeds  in  exactly  the  same  manner,  and  the  same  so¬ 
lution  is  obtained  except  for  the  fact  that  the  Jn(0,t)  term  in 
Eq.  (B13)  is  replaced  by  7n(0,0,  so  that  the  details  are  omit¬ 
ted. 
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